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ABSTRACT 



We present infrared observations of the supernova remnant G2 1.5-0. 9 with the Very Large Telescope, the Canada-France- Hawaii 
Telescope and the Spitzer Space Telescope. Using the VLT/ISAAC camera equipped with a narrow-band [Fe II] 1.64 /jm filter the 
entire pulsar wind nebula in SNR G2 1.5-0. 9 was imaged. This led to detection of iron line-emitting material in the shape of a 
broken ring-like structure following the nebula's edge. The detected emission is limb-brightened. We also detect the compact nebula 
surrounding PSR J1833-1034, both through imaging with the CFHT/AOB-KIR instrument {K' band) and the IRAC camera (all bands) 
and also through polarimetric observations performed with VLT/ISAAC (Ks band). The emission from the compact nebula is highly 
polarised with an average value of the linear polarisation fraction P"!^^ =^ 0.47, and the swing of the electric vector across the nebula 
can be observed. The infrared spectrum of the compact nebula can be described as a power law of index q-ir = 0.7 + 0.3, and suggests 
that the spectrum flattens between the infrared and X-ray bands. 

Key words. Infrared: ISM - ISM: supernova remnants: individual: SNR G21.5-0.9 - pulsars: individual: PSR J1833-1034 



1. Introduction 

In the past years it has been shown that the infrared band 
is an important e nergy range for studying supernova rem- 
nants (SNRs) (e.g.,|01iva et alJU"989l:lKoo et al]|2007l : lLe"e et alj 



120091: lTemim etal.j 120101) and pulsar wind nebulae (PWNe) 
(e.g.. iGraham et al.(ll989HTemim et al.ll2006t ISlane eTa!]l2008 : 
IWilhams et alj|2008l) . especially in the cases where the object is 
heavily obscured by dust in the optical band. Infrared observa- 
tions allow for detecting synchrotron emission fro m relativistic 
partic les perva ding the PWN, e.g. the Crab Nebula (iTemi m et al 
20061) . 3C 58 (ISlane et alj |2008|) . B0540-69.3 ( Willi ams et al 
2008h . and line-emitting material excited either b y a pass- 



The supernova remnant G2 1.5-0. 9 belongs to the class of 
composite SNRs, characterised by a central pulsar-powered 
component, accompanied by an expanding shell of material 
swept up by the supernova blast wav e. Extensive studies of 
G21.5-0.9 carried out in the radio (e.g. Becker & SzvmkowiakI 



ing sh ock, e.g. G 1 1.2-0.3 dKoo et al 2007), 3C 396 ( iLee et alj 
120091) . RCW 103 dOlivaet alJll989l) or by photoionization by 
the synchrotron con t inuum, e.g. filamen ts in the Crab Nebula 
(iGraham et al.lll989l: iTbrnim et al.ll2006l) . This in turn provides 
essential information about the dynamical evolution and physi- 
cal conditions in these objects. Infrared emission may also in- 
clude a component of dust continuum, giving information on 
dust in the remnant v icinity or freshly synthesized in the ejecta 
(IWilUams et al.ll2008h . Moreover, polarisation studies in the op- 
tical or near infrared may be an effective tool for studying the 
inner regions of PWNe as p r oposed by recent MH D simula- 
tions (iBucciantini et alJl2005t iDel Zanna et al.ll2006h . Through 
polarimetry we can learn about the magnetic field structure and 
the flow speed in the wind termination shock area. 
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119811: iFurst et al.|l9 88':'Bietenholz & Barter2008') and in X-ravs 
(e.g. ISlane e t al.| |2000: Safi- Harb et al. 2001; Bocchin oetaD 
120051: iMatheson & Safi-Harbll2010l) reveal the complex struc- 
ture of this object. It consists of a radio and X-ray bright pul- 
sar wind nebula of radius ~ 40" which is suiTounded by a dif- 
fuse and faint X-ray halo, due to dust scattering, with a brighter 
limb of ra dius ~ 150", which traces the shell of the associ- 
ated SNR dBocchino et alJl2005h . A very X-ray bright compact 
source (radius ~ 2") is situated in the centre of the system. It 
harbours one of the most energetic pulsars in our Galaxy with 
a spin -down lu minosity of E = 3.3 x 10^^ erg s ' dGupta et al.l 
,200l ICamilo' et al. 2006i) . The pulsar J1833-1034 has a period 
of 61.8 ms and a characteristic age of 4800 years. A recent age 
estimate based on a measurement of the PW N expansion rate 
in the radio band (iBietenholz & Barteill2008h places this pulsar 
among the youngest in our Galaxy, with an age of 870^j5q yr 
(or even lower, in the case of decelerated expansion), much less 
than the characteristic spin-down time. The distance to the sys- 
tem is estimated, based on HI and '•'CO measurements, to be 4.8 
kpc (Tian & Leahv 2008) with an uncertainty of 0.4 kpc (Tian, 
personal communication). Thus, SNR G2 1.5-0. 9 is one of only 
a few "complete" systems in which a shell supernova remnant 
is observed to contain a pulsar and pulsar-wind nebula. In some 
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composite SNRs, the pulsar is not directly observed; in some 
young PWNe, a shell is either not observed or can be separated 
from the PWN only with difficulty. An object very similar to 
G2 1.5-0. 9 is B0540-69.3, with a similar age and pulsar spin- 
down luminosity, but its location in the LMC (50 kpc distant) 
makes spatial resolution of the PWN much more difficult. 

Due to its position in the inner Galaxy, G2 1.5-0. 9 is heav- 
ily obscured in the optical band, where the estimated value of 
interstellar extinction deriv ed from the hydro gen column den- 
sity A^H - 2 2 x 10 ^^ cm ^ dSlane et alJIiOOOh reaches Ay ^ 10 
dGorensteinl Il975h . The interstellar extin ction drops dramati - 
cally when moving to longer wavelengths dCardelli et al.ll 19891) . 
yielding an estimated near-infrared extinctions A/f ^ 1.2 and 
Afj ^1.9 toward G21.5-0.9. The first mid-infrared observations 
of the remnant taken with the ISOCAM instrument on-board 
the ISO observatory show extended emission w hich matches the 
morphology of the PWN, most clearly at 15 jjm dGallant & Tuffsl 
Il999l) . If interpreted as synchrotron emission, the extracted in- 
frare d fluxes fitted to gether with the previously existing ra- 
dio ( Salter et al. 19891 and references therein) a nd X-ray data 
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1 986t lAsaoka & Kovamal 1 1 9901) by a broken 
power-law model would imply a spectral index of ax - 10 
(Fy cc E^"^) between the infrared and X-ray fluxes. 

SNR G2 1.5-0. 9 was observed in the near-infrared band 
with the Canada-France-Hawaii Telescope, the ESO Very Large 
Telescope and the Spitzer Space Telescope with the aim of de- 
tecting synchrotron emission originating from relativistic par- 
ticles accelerating in the PWN magnetic field, and detecting 
clumped supernova ejecta and/or filaments. A detailed descrip- 
tion of the observations and data analysis can be found in Sect.|2l 
We present results in Sect. |3] their interpretation in Sect. |4]and 
conclusions in Sect. |5] 

2. Observations and data analysis 

2.1. CFHT/AOB-KIR 

Imaging of G2 1.5-0. 9 was carried out using the Canada-France- 
Hawaii Telescope equipped with the Adaptive Optics Bonnette 
(AOB) and the KIR near-infrared camera (together named AOB- 
KIFQ), in order both to look for fine structure in the synchrotron 
nebula and to mitigate confusion with stars in this crowded 
Galactic field. The observations were performed in the K' fil- 
ter, a slightly shorter-wavelength version of the standard K filter 
designed to reduce atmospheric background. The details of the 
obtained data can be found in Table[T] Basic data reduction was 
carried out using the noao.imred.ccdred package of the IRAF0 
Data Reduction and Analysis Facility. For each individual expo- 
sure, bias and dark were subtracted, the flat-field correction was 
applied using dome flats and bad pixels were suppressed. 

Due to the small field of view of the instrument (36" x 36") 
as compared to the expected size of the source, successive obser- 
vations were taken following a raster pattern on the sky so as to 
cover the desired field. In order to reconstruct the full image of 
the sky covered by the observations a mosaicking procedure us- 
ing the IDL programming language was developed. The proce- 
dure makes use of the IDL Astronomy User's Library functions 
CORRELJMAGES and CORRM AT ANALYZE to determine rel- 
ative shifts between subimages with 1 pixel accuracy. A sky level 



' http://www. cfht.hawaii.edu/Instmments/Detectors/IR/KIR/ 
^ IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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Fig. 1. Top panel: CFHT/AOB-KIR mosaicked image of the central 
part of SNR G21. 5-0.9 taken in K' band, smoothed with 2-dimensional 
Gaussian with a FWHM = 0. 1 1" . Bottom panel: the same sky region as 
in the top panel imaged with VLT/ISAAC in Ks band. The image shown 
is the averaged total intensity / (see Sect. l2.2t . For both panels a star de- 
picts the position of PSR J183 3-1034, while the ex tent of the compact 
core seen with Chandra HRC iCamilo et ai]l200d) is shown with con- 
tours. The observations taken with CFHT/AOB-KIR have higher reso- 
lution than the ones obtained with VLT/ISAAC, which has direct im- 
plications for a flux detemiination of an extended emission structure 
visible around the pulsar position (see Sect.|3). 



for each of the subexposures is first computed and subtracted us- 
ing IDL DAOPHOT SKY procedure. (This has the consequence 
that any large-scale emission extending smoothly over most of 
the 36" X 36" subexposure field of view would be subtracted 
out; any such large-scale synchrotron emission could however 
be detectable in polarisation with VLT/ISAAC, see Sect. 12.21 ) 
Having determined the offsets, in the final step a mosaic (size of 
~ 50" X 60") is created by stacking up and averaging the offset 
and sky corrected subimages. The top panel of Fig. [1] shows the 
final mosaicked image. The presented region is smaller than the 
size of the PWN in SN R G21. 5-0.9 and co vers the X-ray com- 
pact nebula region (e.g. lCamiloetal.ll2006l) . 



2 



A. Zajczyk et al.: IR imaging and polarimetric observations of tlie PWN in SNR G2 1.5-0.9 



Table 1. The CFHT/AOB-KIR and VLT/ISAAC observations of 
G21.5-0.9 



Instmment 


Number 


Exposure 


Filter 


Filter 


Date 


name 


of 


time per frame 


central A 


width A/l 


of 




exposures 


[s] 


[^im\ 


[}im\ 


observations 


CFHT 


AOB-KIR 


31 


120.0 


2.115 


0.350 


02.08.2001 


AOB-KIR 


54 


120.0 


2.115 


0.350 


03.08.2001 


VLT 


ISAAC-POL" 












x = 2r 


9 


20.0 


2.160 


0.270 


07.09.2002 




27 


20.0 


2.160 


0.270 


08.09.2002 


X = 1T 


2 


20.0 


2.160 


0.270 


07.09.2002 




34 


20.0 


2.160 


0.270 


08.09.2002 


ISAAC-/MG* 














9 


68.0 


1.64 


0.025 


08.07.2002 




9 


68.0 


1.64 


0.025 


22.07.2002 




10 


68.0 


1.71 


0.026 


08.07.2002 




10 


68.0 


1.71 


0.026 


22.07.2002 



Notes. Polarisation mode of ISAAC instrument (Sect. 12.2) 

Imaging mode of ISAAC instrument (Sect. l2!4] l 



Within a circular region of 15" in diameter centred on the 
pulsar p osition, stars which h ave their brightness in the 2MASS 
K band dSkrutskie et al.ll2006h determined were identified. Their 
fluxes, measured via aperture photometry, allowed for calcula- 
tion of a photometric zero point between the instrumental sys- 
tem and the 2MASS K band. The photometric zero point was 
later used in the transformation of the measured flux to an ap- 
parent magnitude for the extended emission structure detected 
around the pulsar position. 

2.2. VLT/ISAAC - polarimetry 

Polarimetric observations of G2 1.5-0. 9 were taken in the Ks fil- 
ter with the IS AACQ instrument mounted on the ESO Very Large 
Telescope. The instrument uses a Wollaston prism, which splits 
the incoming radiation beam into two perpendicularly polarised 
beams, giving on one CCD plate two images of the same object 
separated by ~ 2 1 " . In order to avoid sources overlapping, a spe- 
cial mask made of alternating opaque (width ~ 24") and trans- 
mitting stripes (width ~ 20") is inserted before the Wollaston 
prism. This instrumental set-up produces a CCD image contain- 
ing altogether six 20" x 150" subimages of the sky - three differ- 
ent stripes of the sky imaged at the same time in two perpendicu- 
lar polarisations. Due to the non-standard mosaic reconstruction 
required, dedicated procedures performing the basic data reduc- 
tion were developed. A sky flat field was obtained for each of the 
six stripes separately from median combination of all the images 
taken with the same angle x (see below), each normalised to its 
median. A raster pattern was applied when performing the po- 
larimetric imaging of G2 1.5-0. 9 in order to obtain continuous 
coverage of a FOV containing the entire PWN. The complete 
image in any of the polarisation beams was then obtained using 
the mosaicking procedure described in Sect. 12.1] adapted to the 
specific character of the ISAAC polarimetric data, including sky 
level determination for each stripe in each image. 

The polarimetric imaging was carried out with two different 
position angle settings of the instrument x (measured North to 

' http://www.eso.org/sci/facilities/ paranal/instruments/isaac/ 



East), namely with;t' - 27° and;^' — 72°, allowing for determi- 
nation of a linear polarisation degree P]^, which in terms of the 
Stokes parameters /, Q and U (see e.g. ldi Serego AUghieril 19971 
for the Stokes parameters definition) can be expressed as: 

Pl- ^Pl + Pl, (1) 

where: 

Q U 
Pq^j and ^'u = y (2) 

and determination of a polarisation angle: 

6 = 0.5 arctan . (3) 

We note that the instrumental set-up is unable to measure any 
circular polarisation which might also be present. The derived 
polarisation angle (Eq. [3) is expressed in the detector reference 
frame. To refer it to the celestial reference frame a counter- 
rotation of 27° with respect to the instrument position angle 
is app Ued following e.g. Eq. (9) of iLandi DeglTnnocenti et aU 
(l2007h . In the celestial reference frame the polarisation angle is 
measured with respect to the North Celestial Meridian (North 
to East). The bottom panel of Fig. [T] shows the central part of 
G2 1.5-0. 9, where the X-ray compact nebula is situated, as seen 
in the total intensity /. 

Polarisation standard stars were also observed and their lin- 
ear polarisation degree and angle were measured in order to 
verify determination of polarisation parameters for the compact 
nebula region. The standard stars were chosen following the 
ISAAC list of polarised standard stars. We used R CrA15 and 
HD 188112. The measurements for R CrA15 in the Ks filter 
yield Pl = 0.017 + 0.004 and 6» = 95° + 13° determined North 
through East. The result is in fairly good agreement with the 
known polarisation for this object at shorter wavelengths - e.g., 
in the H band Pl.h = 0.0144 + 0.0006 and Gh = 100° + 2° 
and in the J band P , , = 0.0228 + 0.0003 and 6»y = 96° + 1° 
(IWhittet et al.lI992l) . For HD 188 1 12 we get Pl = 0.007 + 0.005 
ande^ 118° + 14°. 
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For the central region of G2 1.5-0. 9 two maps, of linear 
polarisation fraction f l and linear polarisation intensity L - 
yjQ- + U~, were created. They are presented in Figure For 
presentation purposes pixels with signal of Icr below the sky 
value - as determined in the total intensity image / - were masked 
in the Pl map. In both maps an extended blob of polarised emis- 
sion surrounding the PSR J1833-1034 position is visible. In or- 
der to show the behaviour of the polarisation angle across the 
nebula the electric field vectors are overlaid on the L map. They 
were determined from the signal integrated over circular regions 
(r ~ 0.45") spread across the polarised emission nebula. A grad- 
ual change in the polarisation angle is visible when moving from 
the SE to NW part of the torus region and its implications are 
discussed in more detail in Sect. 14.11 

Similarly as for the CFHT observations, the flux of the emis- 
sion coincident with the putative pulsar wind torus identified in 
the averaged total intensity image / was determined. To compute 
the photometric zero point, flux measurements for the polari- 
metric standards were used together with their 2MASS K magni- 
tudes. 

2.3. Spitzer 

SNR G2 1.5-0. 9 was observed with the Infrared Array Camera 
(IRAC) and Multiband Imaging Photometer (MIPS) onboard 
the Spitzer Space Telescope as part of observing programme 
ID 3647. Using the Spitzer Heritage Archive, po^r-Basic 
Calibrated Data (post-BCDs) containing mosaicked images of 
the G2 1.5-0. 9 region were obtained for all IRAC bands and the 
24 jum MIPS band. The IRAC observations at 3.6, 4.5, 5.8 and 
8.0 fim were carried out on 15 September 2005, while the MIPS 
24 fim observation was taken on 1 1 April 2005. Figure |3] shows 
flie central part of G2 1.5-0.9 in all IRAC bands. 

2.4. VLT/ISAAC - [Fe II] imaging 

G2 1.5-0. 9 was also observed in the [Fe II] 1.64 yum filter and in a 
neighbouring narrow -band 1.71 //m filter with the ISAAC instru- 
ment (in standard imaging mode). A dark correction was applied 
to all frames following the standard approach using IDL proce- 
dures written to handle the ISAAC data. In each of the filters the 
successive exposures were performed in a wide raster pattern so 
as to cover the sky area where emission from the PWN, as well 
as a dense part of the SNR shell, could be expected. Large shifts 
between exposures were used to ensure that a given camera pixel 
fell outside any diffuse emission from G2 1.5-0. 9 in most frames. 
A sky flat field for each of the narrow-band filters could then 
be constructed by median combination of all exposures taken 
on the same night (each scaled such that its median was unity). 
These sky flat fields were divided into the science frames follow- 
ing the standard approach. The IDL mosaicking procedure that 
also determines the sky level for each of the subexposures (see 
Sect. 12. II ) was applied to obtain the final image in both the [Fe II] 
1.64 yum and 1.71 pm filters. 

In order to minimise confusion from faint stars which is evi- 
dent in the [Fe II] 1 .64 pm image (Fig.|4] panel a), we subtracted 
from it the 1.71 pm image (Fig. H] panel b). To compensate 
for the difference in stellar fluxes in different filters, an overall 
normalisation factor was computed from measured star fluxes 
in both nan^ow-band filters (/^i.64/jm/^'i.7i//m ~ 0.84) assuming 
that star colour effects can be neglected over the small wave- 
length difference. A simple subtraction of the aligned and scaled 
1.71 pm image yields rather conspicuous stellar residua, how- 



a) 




' L 

I 



[ 



35.7^ ]8^53"'55.5^ 33.3' 

Right oscension 



b) 




L 

[ 



35..-* 18^33^33,5^ 53.3^ 

Right oscension 



Fig. 2. a) Linear polarisation degree Pl map of the compact core re- 
gion of SNR G2 1.5-0.9. The map was constructed using polarimetric 
imaging observations performed with VLT/ISAAC and was smoothed 
using a 2-dimensional Gaussian with FWHM = 0.45". Colour-coded 
is the measured value of Pj^. Con tours show the exte nt of the com- 
pact core seen with Chandra HRC iCamilo et alj|2003) . b) Map of lin- 
early polarised intensity L = -^/Q^ + smoothed with a 2-dimensional 
Gaussian with a FWHM = 0.45". Colour-coded is the measured value 
of L. White bars represent electric field vectors with their length propor- 
tional to the linear polarisation fraction /'l and their orientation defined 
by the polarisation angle 9. A vertical bar placed in the bottom right 
comer, for reference, shows an electric field vector with Pl = 1 and 
9 = 0°. In both panels a star shows the pulsar's position. 



ever, due in large part to different PSFs in the two images owing 
to changing observing conditions. An improved subtraction was 
thus performed using a PSF equalisation procedure, which also 
ensured image alignment to sub-pixel accuracy, proceeding in 
the following way. We used 20 copies of the two original im- 
ages with integer pixel offsets (all those with a total 2-D offset 
less than 2.5 pixels), and computed the linear combination of all 
these images, and a uniform background, which yielded minimal 
residua and a minimal contribution from the offset images. This 
optimisation was performed using a central image area which 
had negligible diffuse emission but a PSF as close as possible 
to the average one within the PWN area. The result is shown in 
panel c of Fig. H) the effectiveness of the subtraction procedure, 
in the PWN area in particular, can be gauged by comparison 
with panel a, which shows the unsubtracted image in the same 
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Fig. 3. IRAC images of the central region of SNR G21.5-0.9: 3.6 fim 
- top left, 4.5 yum - top right, 5.8 yum - bottom left and 8.0 /jm - bottom 
right. All panels cover the same field of view. Contours in the 4.5 yum 
image show the extent o f the compact nebula as seen with Chandra 
HRC iCamilo et al1l2006h . The star marks the position of PSR J1833- 
1034. In the top left panel circles depict: large central circle - flux ex- 
traction region for the compact nebula; smaller circles: regions where 
the background level was estimated. These regions were used for mea- 
suring the flux in all IRAC bands (see Sect. 13. 2t . Each image was 
smoothed with a 2-dimensional Gaussian with FWHM = 1.8". 



intensity scale. Some residua remain, including "halos" around 
brighter stars, and may be due in part to non-linearities in in- 
strument response, as well as to intrinsic differences in stellar 
spectra between the 1.64 //m and 1.71 /im bands. 

Panel c of Figure |4]reveals, after the above subtraction, dif- 
fuse emission in the [Fe II] 1 .64 fim filter coincident with parts 
of the PWN in G2 1.5-0. 9. Extended emission in the shape of a 
partial ring is clearly visible. Chandra ACIS contours are over- 
laid to show the PWN extent in X-rays. The surface brightness 
for the NE and S rims of the [Fe II] 1 .64 fim emission was com- 
puted. The flux calibration was performed using stars that have 
2MASS H magnitudes determined and that could be identified 
in the observed field. Detailed discussion of the results can be 
found in Sect. 13. 31 

3. Results 

3.1. Compact Core - near-infrared 

The field of view {FOV) of the central part of the AOB-KIR mo- 
saicked image with the best S/N is approximately 36" x 36". 
It covers the central part of the PWN as identified by X-ray 
observations, and i n part icular the compact nebula around the 
pulsar (jsiane et al] |2000|) . The observed field is rich in stars 
which can confuse any NIR emission coming from G21.5- 
0.9. Using a brightness scale in which stars are saturated, a 
clearly extended region of emission (compared to the instru- 
ment's FW//Maob-kir - 0.24") emerges from the background 
(Fig. [T] top panel). Its size of ~ 4" in diameter is comparable 
with the size of the X-ray compact nebula. Its apparent centre 
is slightly shifted towards the SW with respect to the pulsar po- 



sition a = 18'^33'"33':57 and 6 = -10°34'07'.'5 (ICamilo et al.l 
l2006h . but more extended emission may be confused with the 
stars that are situated to the NE of the pulsar It is worth not- 
ing that the X-ray morphology of the emission region around 
t he p ulsar position jCamilo et al.ll2006t iMatheson & Safi-Harbl 
2010h is similar to the one seen in the present near-infrared data, 
and can be caused by a Doppler-boosted toroidal structure sur- 
rounding the pulsar 

To extract the flux from the extended emission nebula a cir- 
cular aperture 6" in diameter around the pulsar position was 
used. The contributions to the flux of the nebula from stars lo- 
cated within a 5" radius of the pulsar were subtracted by per- 
forming aperture photometry. The measured flux for the com- 
pact nebula corrected for interstellar extinction is Fk' - (0.82 + 
0.60) mJy. The existence of this emission blob is confirmed by 
the polarimetric observations. No other significant diffuse emis- 
sion is detected. 

Each of the mosaicked images, i(0°}, 1(45°), i(90°) and 
i(135°), obtained from the subimages taken in the polarimetric 
mode of the ISAAC instrument, has a FOV of ~ 100" x 200". 
In all cases it includes the entire region where the X-ray PWN 
is identified. The detected diffuse emission, which can be seen 
in the total intensity image / (Fig.lT] bottom panel), is restricted 
only to the very central region of the PWN. It is slightly extended 
and compatible in size and position with the emission blob de- 
tected with AOB-KIR in the K' filter 

Measurements of the degree of linear polarisation f l and the 
polarisation angle 9 reveal highly polarised emission from the 
nebula with a regular structure of the electric field (Fig. panel 
b). As can be seen from the Pl map (Fig. |2l panel a) this po- 
larised emission extends within the X-ray compact nebula and 
is only important where the star contamination is minimal. To 
the NE of the pulsar position, but still within the X-ray compact 
nebula, a region where the linear polarisation degree is smaller 
than 0.1 can be identified (Fig. |2l panel a). As seen from the 
total intensity image / (Fig.[T] bottom panel), this is also the re- 
gion where strong stellar contamination is present. The linear 
polarisation fraction averaged over the emission region is high, 
pfl>,? ^ Q 4-7 _j_ 0.23; moreover, in some regions within the polari- 
sation nebula Pl is even as high as 0.6. The high degree of linear 
polarisation points to the synchrotron nature of the observed ra- 
diation and a highly ordered magnetic field structure. 

As mentioned before, the measurements of 6 across the com- 
pact nebula reveal a regular structure of the radiation electric 
field which is illustrated as white vectors in the b panel of Fig.|2] 
In the central part of the nebula the electric field vectors are ori- 
ented in the NE-SW direction. In the outer regions of the neb- 
ula, vectors tend to bend away from the NE-SW direction in a 
twofold way: moving to the NW vectors rotate anticlockwise; in 
the SE part vectors rotate clockwise with respect to the electric 
field direction determined in the central part of the polarisation 
nebula. The change in orientation of the electric field vectors is 
accompanied by a steady decrease in the linear polarisation frac- 
tion Pl- 

A circular aperture similar in size to that used for the CFHT 
flux measurements was used to extract the counts from the de- 
tected polarised emission blob. Star contamination to the mea- 
sured flux was reduced by their standard photometry measure- 
ments. The dereddened flux for the emission blob derived from 
the averaged total intensity is F^s - (1.28 + 0.84) mJy; although 
it is compatible within uncertainties, we consider this value less 
reliable than the one obtained from the AOB-KIR observations 
(see Sect. 14211. 
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Fig. 4. a) ISAAC mosaicked image of SNR G2 1.5-0.9 in the [Fe II] 1.64 yum naiTow-band filter. The full field of view is presented to show also 
the North Spur region, b) Same field of view as in panel a, but imaged in the 1.71 /jm narrow-band filter. Faint diffuse emission at the NE and SW 
of the PWN is visible in 1.64 jim filter, and not in 1.71 yum. c) The [Fe II] 1.64 yum image after subtraction of the image in the 1.71 yum narrow-band 
filter. Some stellar residua are still visible; the net flux in these residua is typically around 9% (RMS) of the coiTesponding stellar flux. For all 
panel s: the star depicts the position of PSR J1833-1034; contours represent the extent of the X-ray PWN seen by Chandra ACIS jSlane et al] 
l20Q(ih . All images are presented with the same intensity scale. For presentation purposes, images were smoothed with a 2-dimensional Gaussian 
with cr = 0.28". The images shown include the data from the first night of observation (08.07.2002) only, in order to yield more uniform coverage 
of the PWN and North Spur regions. 
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Table 2. Spitzer IRAC and MIPS observed and extinction 
corrected fluxes of the compact nebula. Correction for inter- 
ste Uar extinction was per formed using the local ISM table 
of |Ch lar & TielensI (l2006l) . where Ak - 1-15 wa s estim ated 
from A^H = 2 2 x l O^^cm'^ following fGorensteinl (1 19751) and 
ICardelUet al.]([T98l . 



Instrument 


A 
[jum] 


Observed Flux 
[miy] 


De-reddened Flux 
[mJy] 


IRAC 










3.6 


0.53+0.19 


0.94 ± 0.45 




4.5 


0.74 ± 0.14 


1.19 ±0.49 




5.8 


0.88 ± 0.20 


1.34 ± 0.49 




8.0 


0.89 ± 0.22 


1.38 + 0.71 


MIPS 










24.0 


2.47 ± 0.67 


4.37 + 2.13 



sion associated with the PWN. Its shape follows the outer con- 
tours of the X-ray nebula (Fig. |4] panel b). The brightness dis- 
tribution is non-uniform; a partial ring-like structure has a max- 
imum of emission in the NE part. Going towards the SW, the 
emission fades leaving only a few faint patches of [Fe II] emis- 
sion which fall within the X-ray contours of the PWN (Fig. |4] 
panel c). The brightest NE portion of the ring shows clear limb- 
brightening with a well-defined outer edge. To the SW from the 
pulsar position an emission blob, probably an artefact associ- 
ated with the bright star (as seen next to other bright stars in 
the FOV) rather than with the PWN material, is apparent. The 
mean surface brightness of the NE part of [Fe II] line emission 
is 5[Fe ii]i.64/im = (3.3 ± 0.5) X 10"'* erg s"' cm"^ arcsec"- while 
for the southern rim it is ^[Fe u]i.64/jm - (1-7 + 0.3) x 10"'* erg 
s"' cm"- arcsec"^. 



3.2. Compact Core - mid-infrared 

No large-scale extended emission associated with the PWN is 
seen in IRAC bands. However, in the PWN centre a compact 
emission blob is present (Fig.|3]l. It is clearly visible in all IRAC 
bands and its extent is comparab le with the size of t he putative 
pulsar wind torus seen in X-rays (ICamilo et al.ll2006l) and in the 
near-infrared. For the IRAC observations the source flux was ob- 
tained using an extraction region with a radius of 3 .6" (large cir- 
cle in the top left panel of Fig. 13). To estimate the background 
level, circular regions with a 2" radius (small circles in the top 
left panel of Fig. |3]l were chosen in the vicinity of the source ex- 
traction region. The average value from these regions was taken 
to represent the background level at the source position. Because 
the aperture used for the flux estimation in the IRAC bands was 
different from the aperture used for the Spitzer calibration stars 
an aperture correction, as described in IRAC HandbookQ, was 
apphed to the computed fluxes. In order to estimate the compact 
nebula's flux at MIPS 24 pm standard aperture photometry was 
used. A source extraction region with radius of 7" was used. The 
background was estimated using an annulus with an inner radius 
of 7" and an outer radius of 13". To obtain the final value of the 
measured flux an aperture correction was applied as described in 
the MIPS HandboolQ. Spitzer IRAC and MIPS 24 pm fluxes are 
presented in Table |2] The flux uncertainties reflect the variations 
in the background level between different background regions, 
as well as the uncertainty in the extinction correction for the de- 
reddened fluxes. 



3.3. Pulsar Wind Nebula 

The mosaicked image taken in the [Fe II] 1.64 pm and 1.71 pm 
filters (see Sect. covers a FOV of ~ 230" x 160". This 
includes both the PWN region and the North Spur region (a 
bright spot located in the northern part o f the X-ray halo sur- 
rounding the PWN (iBocchino et al.ll2005h ). By subtracting the 
1 .7 1 pm image (corrected for the PSF and average star flux dif- 
ference between the filters) from the [Fe II] 1 .64 pm image (see 
Sect. 12.4b we minimise the stellar contamination to any [Fe II] 
1.64 pm line emission that can be present in the observed re- 
gion. The resultant image is presented in panel c of Fig. |4] No 
[Fe II] line emission associated with the North Spur is detected. 
There is, however, extended (~ 80" in diameter) [Fe II] emis- 



4. Interpretation 

Both the emission detected with AOB-KIR, Spitzer and 
the ISAAC polarimetric observations and also the structure 
seen in the [Fe II] 1.64 pm emission line, can be in- 
terpreted within a composite supernova remnant evolution- 
ary scheme (e.g. Gaensler & Slanell2006t Ivan der Swaluw et al.l 
l200lllRevnolds & Chevaliedl984l) . According to this evolution- 
ary scheme a supernova (SN) explosion drives a blast wave into 
the surrounding interstellar medium (ISM). At the same time in 
the interior of the expanding SNR a pulsar wind nebula evolves. 
It is powered by the newly bom pulsar through its relativis- 
tic magneto-hydrodynamic wind. At early times, the pulsar en- 
ergy input causes the PWN's expansion to accelerate, driving 
a shock into the inner edge of the uniformly expanding ejecta. 
Inside the PWN, at the point where the pulsar wind ram pres- 
sure is equilibrated by the pressure exerted by the surrounding 
medium, a pulsar wind termination shock (WTS) is formed. The 
schematic structure of such a young composite SNR is presented 
in Figure [5] 



4.1. Polarisation - Wind Termination Shock 

The high degree of linear polarisation measured across the emis- 
sion blob seen in the Ks band (see Sect. 13.1b points to its syn- 
chrotron origin. Simultaneously, the orientation of the electric 
field vectors (Fig. |2] bottom panel) suggests that the region 
where synchrotron radiation is produced has a toroidal magnetic 
field. Such conditions are expected to be found in the wind ter- 
mination shock region (Fig. |5]l- The relativistic particles flow- 
ing upstream of the shock do not radiate. Only after being ac- 
celerated at the shock, can they produce synchrotron radiation 
in the downstream flow. The best exam ple in suppor t of such a 
scenario is the torus in the Crab Nebula (lHesteril2008[ and refer- 
ences therein). When viewed in polarised light it shows magnetic 
field structure that is aligned along wisps and fibrous structures 
which are present in the torus, and which follow its shape. For 
some of the wisps a high degree of Unear polarisation, near 0.7, 
is measured. 

Recent MHD simulations studying optical polarisation in 



http://ssc.spitzer.caltech.edu/irac/iracinstrumenthandbook/ 
^ http://ssc.spitzer.caltech.edu/mips/mipsinstrumenthandbook/ 



the inner parts of pulsar wind nebiilae ( Bucciantini et al.l 120051: 
iDel Zanna et al J 120061: IVolpi et alJl2009h predict that emission 
from a wind termination shock region should be highly po- 
larised. In these simulations a pure toroidal structure of magnetic 
field and a radial flow of paiticles are assumed. The maximum 
value of linear polarisation fraction, close to 70%, is expected 
from the WTS region close to the torus symmetry axis as pro- 
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Fig. 5. Schiematic view of a young composite supernova remnant struc- 
ture (not to scale). The supernova (SN) blast wave propagates in the 
interstellar medium (ISM). Inside the SN remnant a pulsar wind neb- 
ula (PWN) powered by the pulsar (star) is formed. The PWN ex- 
pands into the cold SN ejecta, which from the outer side are being 
heated/compressed by the reverse shock facing inwards in the SN rem- 
nant. Inside the PWN, where the pulsar wind pressure balances the 
pressure of the suiTounding matter a wind termination shock (WTS) 
is formed. The inner solid line indicates the shock driven into the in- 
ner ejecta by the accelerating PWN. Dotted line depicts the position of 
the contact discontinuity separating the shocked SN ejecta overtaken by 
the PWN outer shock, and the shocked pulsar wind plasma having gone 
through the WTS. Detailed discussion of SNR/PWN structure can be 
found in e.g. lKennel & Coronitil ( fT984h or lTruelove & McKed < 19991) . 



jected on the sky. When moving away from the axis, the polar- 
isation progressively decreases to values smaller than 50%. In 
addition, close to the torus edges, regions tha t are almost com- 
pletel y depolarised can be found (see Fig. 2 of lBucciantini et aP 
l2005h . Simulations also show that the polarisation angle should 
change across the emitting region. In the parts close to the sym- 
metry axis of the torus, the electric vector should be aligned with 
the axis and its inclination with respect to the axis should in- 
crease while moving away from it. The region where the electric 
vector becomes perpendicular to the symmetry axis depends on 
the inclination angle of the torus with respect to the sky, and also 
on the bulk flow speed. Concerning the latter dependence, the 
position where the polarisation vector becomes perpendicular to 
the axis moves closer towards the axis with increasing value of 
the bulk speed (iBucciantini et al.ll2005h . 

Similar behaviour to that predicted by MHD simulations 
is found in our polarisation observations of the central part of 
G2 1.5-0. 9 (see Sect. 13. Il l, which supports the claim that in the 
region of the wind termination shock a toroidal magnetic field is 
present. From the investigation of the observed polarisation an- 
gle swing, the probable orientation of the symmetry axis of the 
toroidal field as projected on the sky is *P ~ 40° (measured North 
to East). In this scenario, is also the position angle of the pul- 
sar spin axis as projec ted on the sky, whic h is in agreement to 
the value proposed bv lCamilo et al.l (l2006h . The nebula seen in 
the Ks band can be described as an ellipse with semimajor axis 
of 01 ^ 2.9" and semiminor axis of 62 - 1.8". The ellipse mi- 
nor axis is oriented at ^', so it is parallel to the toroidal magnetic 
field symmetry axis as projected on the sky. However, due to the 
Doppler boosting experienced by the relativistic particles escap- 
ing the wind termination shock, the visible size of the emission 
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Fig. 6. Spectral energy distribution of the compact nebula. The dia- 
mond shows the flux obtained with CFHT/AOB-KIR. Spitzer measure- 
ments obtained in the IRAC 3.6, 4.5, 5.6, 8.0 yum bands and the MIPS 
24 yum band are presented with triangles. All the infrared fluxes are cor- 
rected for interstellar extinction. Light grey bow tie region shows the 
uncertainty of the power-law fit to the infrared fluxes. The fitted IR pho- 
ton index is Fir = 1.7 ± 0.3. The black bow tie represents the intrinsic 
X-ray spectrum (with its unc ertainty) of the comp act nebula as obtained 
from Chandra observations dCamilo et al.ll2006n . It is described by the 
photon index T = 1.5 ± 0.1 (see text). 



region cannot be simply translated into the shock radius. A vital 
part of the WTS may be invisible because the size of the emit- 
ting re gion highly depends on the bulk flow speed in the WTS 
region (IBucciantini et al ] l2005h . If the emitting torus is not seen 
edge-on, but is inclined with respect to the plane of the sky, the 
Doppler boosting effect could possibly explain the shift of the 
compact nebula emission towards SW with respect to the pul- 
sar position. To more accurately estimate the orientation of the 
toroidal magnetic field and the flow speed, detailed modelling 
is needed . However, by siinple co r nparison with the nume rical 
results of Buccian tini et al.l (120051) . iDel Zanna et al.l (|2006|) and 
Vol pi et al.. (2009) one may expect that the particle flow speed in 
the WTS region is high (> 0.2c). 

4.2. Nonthermal Emission Spectrum 

Figure |6] presents the spectral energy distribution of the com- 
pact nebula based on our flux measurements in the CFHT/AOB- 
KIR, Spitzer IRAC 3.6, 4.5, 5.6, 8.0 fim and MIPS 24 fim bands, 
and its intrinsic X -ray spectrum obtain ed from Chandra observa- 
tions presented bv lCamilo et al.l(l2006l) . The plotted IR fluxes are 
the values cor rected fo r the interstellar extinction using the local 
ISM table of iChiar & T ielens (2006), where Ak ^ 1.15 was es- 
tim ated from A^H 7 2.2 x lO^^cm"- following lOorensteinl (1 19751) 
and lCardeniet al.1 (119891) . A simple power law Fir - AE 
where aiR is the spectral index, was fitted to the infrared fluxes 
using a least-squares fitting method that takes into account the 
uncertainties of the measured fluxes as their weights. In the fit 
the flux value obtained from the VLT/ISAAC-POL averaged to- 
tal intensity / was not included because the flux value from 
CFHT/AOB-KIR observations is a more accurate representa- 
tion of the compact nebula flux at ~ 2.1/im. This is because 
the higher resolution (smaller PSF) of CFHT/AOB-KIR allows 
for much more accurate stellar subtraction in the crowded com- 
pact nebula field, which is difficult in the case of VLT/ISAAC- 
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POL and in turn leads to an overestimation of the measured flux. 
The fitting to the rest of the IR fluxes yields the spectral index 
aiR - 0.7 + 0.3. The obtained power law with its uncertainty is 
presented in Fig.|6]as the light grey bow tie region. 

The X-ray spectrum in the 0.5 - 10 ke V energy range was 
obtained from Chandra ACIS observations jCamilo et al.l 20061) 
using a circular extraction region with 4" radius around the PSR. 
The background level was estimated using an annulus with ra- 
dius 4"- 10". Usi ng the best-fit valu e of A^h = 2.2x10" cm"^ for 
the entire nebula dSlane et al.ll2000h this yields a photon index of 
r = 1.46 + 0.02 and unabsorbed flux in the 2.5 - 10 keV energy 
range /^2.5-iokeV = (7.8 +0.2) x 10"'^ erg s"' cm"^. Allowing the 
column density to be determined just from the fit of the compact 
region yields similar values of A^h = (2.3 + 0.1) x 10^^ cm 



F = 



7.7 



+0.5 



X 10 erg s ' cm ^. 



1.56 + 0.04 and f 2.5-iokeV = ' ■ ' _o.i 
Comparison of the fitted IR spectrum with the X-rays (Fig. |6) 
suggests that the spectrum must flatten between infrared and 
X-rays. Given the unusual implications of such a result on the 
underlying particle spectrum, it is important to consider possi- 
ble factors that could complicate the interpretation of the X-ray 
spectrum. There are two primary sources of potential bias: con- 
tamination from the central pulsar, and dust scattering effects. 

Regarding the pulsar contamination, to estimate a point 
sou rce contribution to th e unabsorbed flux of the compact neb- 
ula, ICamilo et al.l (l2006l) use Chandra HRC data to model the 
spatial emission around PSR J1833-1034 as a point source and 
a 2-D Gaussian. Converting their best-fit count rate of the point- 
like component into a flux using Nh = 2.2 x 10^^ cm"- and 
an assumed photon index of 1.5 yields an unabsorbed flux 



1.8 X 10"' erg s"' cm"^ (about 2% of the unab- 



rPSR 

2.5-iokeV 

sorbed flux obtained for the compact nebula). This contribution 
is insufficient to account for the apparent difference in spectral 
index between the X-ray and IR bands (and the reduced flux 
upon subtracting the point source contribution marginally in- 
creases the difficulty in matching the spectra). 

The properti es of X-ray scattering from dust grains are well 
known (see e.g. I Smith & Dwekl[l998l) . but modelling of X-ray 
halos of diffuse sources is a difficult task, because the effect is 
non-local. Bocchino et al. (2005) have modelled the X-ray halo 
detected around G2 1.5-0. 9, showing that it is mostly an effect 
of dust scattering, deriving halo parameters, like the scattering 
optical depth Tjcatt - from this energy dependence it 

is apparent that lower-energy photons are scattered with much 
higher efficiency. An additional effect of scattering is a radial 
dependence of the best-fit column density (A^h), which in G21.5- 
0.9 ranges from ~ 1.5 x 10"^ cm"^ in the outskirts of the PWN, 

up to 

iBocchino et"ai]|2005h 

A consequence of scattering is a deficiency of soft photons in 
the central regions, which leads standard spectral analysis rou- 
tines to overestimate the absorption column density; conversely, 
the surrounding regions receive the scattered photons and show 
an excess of soft photons that leads to a negative bias on the 
estimated A^h- Since in principle similar biases could affect the 
estimates of the spectral index and unabsorbed flux, depending 
on the details of the data analysis, we set upper limits to those 
possible biases. We have chosen the energy range 2.5 - 10 keV 
for our fits to the spectrum of the inner torus: this is a good 
compromise for weighting harder photons (only weakly affected 
by scattering) but still retaining good statistics. With this energy 
cut we cannot independently determine the value of A^h, but in- 
stead consider a range of A^h parameters, chosen from 2.0 to 
2.4 X 10^^ cm"^ (2.2 x 10^^ cm"^ being the average value ob- 
tained from a fit to the whole PWN). 



2.5 X 10^~ cm in the very central regions (see Fig. 3 in 



Table 3. Grid of hydrogen column densities A^h used for deter- 
mination of the X-ray spectrum and the unabsorbed flux in the 
2.5 - 10 keV range. Fits were performed using an annulus with 
radius of 4" - 10" for background level estimation {top part of 
the table), and not using any background region {bottom part of 
the table). 



Nh 
[10^- 



-f2.5-10kcV 

[10"'^ erg s"' cm" 



Background from surrounding annulus 



2.0 
2.2 
2.4 



1 AT, + 004 

-'•^•^ - 0.05 

1 47 + 0.05 

' - 0.04 

1 51 + 

-'•-^^ -0.05 



77 +0.2 
' ■ ' - 0.3 

7 6 + "-^ 

- 0.2 



7.7 



- 0.2 

- 0.2 



No background 



2.0 


1.49 


+ 0.04 
- 0.04 


9.6 


+ 0.2 
- 0.2 


2.2 


1.53 


+ 0.04 
- 0.03 


9.7 


+ 0.2 
- 0.2 


2.4 


1.57 


+ 0.04 
- 0.03 


9.7 


+ 0.2 
- 0.2 



Another question is whether or not to use as background the 
emission of the surrounding PWN: the emission in the surround- 
ings of the torus is partly contaminated by the halo of the torus 
itself, and some of the photons emitted by the outer PWN are 
scattered into the area on which the torus extends. In Table |3] we 
present the results for the grid of cases mentioned above. From 
this we can conclude that conservative error intervals on the pa- 
rameters, taking also into account dust scattering effects, are the 
compact nebula photon index of F = 1 .5 +0. 1 and the unabsorbed 
flux of F2.5-iokeV = (8.7 + 1 .0) X 10"'" erg s"' cm"^. 

We note that a joint fit between the X-ray and IR data, leav- 
ing all parameters free, yields Nh - (2.6 + 0.2) x 10^^ cm"^ and 
F = 1 -76^^ Applying these values to the X-ray data alone pro- 
vides a good fit {x^ - 1.1), though the X-ray fit obtained with Nh 
and F as free parameters is statistically much better {xj = 1.01). 
The results are similar when a dust-scattering component is in- 
cluded in the fitting process. While the spectral modeling thus 
argues for a spectral flattening between the IR and X-ray bands, 
it is clearly important to search for examples of such behavior 
in other sources given that the results here rely exclusively on 
somewhat modest differences in parameters from spectral fitting. 

Even though the torus in G2 1.5-0. 9 is not detected in the 
radio, it would appear that a further break or steepening must 
exist betw een the radio and IR ban ds. Using the radio map 
at 5 GHz dBietenholz & Bartelll2008l) we can estimate an up- 
per limit F^QYi^ - 52 mJy to the radio flux of the compact 
nebula by integrating the signal within a radius of 3" around 
the pulsar position. This can be compared to the 5 GHz flux 
^^ext - 11 Jy which one gets by extrapolating the fitted IR spec- 
trum to the radio band. This extrapolated value Fgxt is much 
higher than the estimated upper limit F^q^^-^, which argues in 
favour of the existence of spectral break between the radio and 
IR bands. This would be expected because the radio spectrum 
of the torus presumably has the same index as the entire PWN 
(whose radio index is flatter than what has been measured here 
for the torus in the IR), because the particles responsible for 
the radio emission will not have suffered significant radiative 
losses in filling the nebula. Similar spectral behaviour is also 



9 



A. Zajczyk et al.: IR imaging and polarimetric observations of tlie PWN in SNR G2 1.5-0.9 



observed for the toras in 3C 58 (ISlane et al 1 l2008h . These two 
cases show that the injection spectrum is complex, and not a 
pure power law. Fleishman & Bietenholz (2007) have calculated 
synchrotron spectra of a power-law distribution of electrons radi- 
ating in magnetic turbulence of prescribed properties. For appro- 
priate choices of parameters, it is possible to produce a spectrum 
that could continue to steepen from radio to X-ray energies, as 
our data imply; however, the suggestion in our data of a flatten- 
ing of the spectrum between IR and X-rays cannot be produced 
by such a model. Moreover, the iFleishman & Bietenholj (l2007h 
calculations were applied to the integrated spectra of the PWNe 
and as such they may not be applicable to explain the spectra 
produced in the vicinity of the wind termination shock, which is 
the case for the spectrum presented in Fig. |6] 

4.3. [Fe II] filaments - PWN shock 

Studies carried out in the past 20 years on infrared emission from 
shell supernova remnants show that [Fe 11] emission seen in shell 
SNRs is either induced by the reverse shock propagating back 
through the fast moving ejecta or produced at the forward shock 
from iron present in the circumstellar and/or interstellar medium. 
In the case of pulsar wind nebulae, by contrast, iron emission, 
if collisionally excited, is produced by the rapidly expanding 
PWN encountering the slow-moving cold ejecta (Fig. |5]l. This 
interpretation has been applied to observed [Fe 11] lines at 17.9 
and 26 urn from G54. 1 +0.3 (Temim et al. 2010) andB0540-69.3 
(IWilhams et al.ll2008h . Our imaging observations demonstrate 
that in SNR G2 1.5-0. 9, the iron emission does in fact originate 
at the location of the PWN-ejecta interaction. The nature of this 
emission then characterizes the otherwise unobservable inner- 
most supernova ejecta, material formed just above the mass cut 
within which material fell back to become the neutron star 

In G2 1.5-0. 9 the detected [Fe II] 1.64 fim emission - hav- 
ing a ring-like structure that is significantly brighter in the NE - 
follows the outer edge of the PWN (Fig.HJi. This behaviour sug- 
gests that the [Fe II] emission is caused by the PWN driven shock 
that compresses the surrounding cold supernova ejecta (Fig.|5]). 
Collisional excitation would then take place in the medium just 
behind the shock front. The low excitation energy means that 
only low shock velocities are required, perhaps as a result of 
shocks being driven into dense clumps in the ejecta, but it is also 
possible that the PWN forward-shock speed is slow due to its 
expansion into moving ejecta. Figuring out which one is more 
plausible in the case of G2 1.5-0. 9 would require a dynamical 
model and/or spectroscopic study in search of density diagnos- 
tics in the line-emitting regions. The observed brightness dif- 
ference between regions in the NE and SW could result from a 
density difference in the medium that the expanding PWN bub- 
ble encounters, which in turn would imply a difference in the 
PWN shock velocity in the NE part with respect to the SW. 

A good example of the PWN expanding into a medium 
of nonuniform density is the Crab Nebula. Such conclusion is 
drawn based on observations of [O III] and [Ne V] emission, 
excited by the PWN shock, which forms so-called skin at the 
boundary of the Crab Nebula. The skin is visible around most of 
the outer edge of the ne bula with the exception of its NW part 
dSankrit & Hestedll997h . The lack of [O III] and [Ne V] emis- 
sion is explained as due to variations in the density of the ejecta 
into which the nebula is expanding. This view is also supported 
by the morphology of Rayleigh-Taylor fingers (LoILet_al..., .2007). 
which shows SE-NW asymmetry. In the NW part of the Crab 
Nebula there are fewer fingers comparing to the SE. This asym- 



metry results f rom a pre sho ck ejecta density lower in the NW 
than in the SE (lHesteijl2008l and references therein). 

However, another possible sc enario may explain th e [Fe II] 
1.64 iim line emission in PWNe. iGraham et al.l (Il989l) showed 
that the near-infrared iron emission detected in the Crab Nebula 
originates from optically thick filaments which are photoion- 
ized by the UV-X-ray component of the synchrotron continuum 
present in the PWN. In the general ISM, a measured line ra- 
tio ([Fe II] 1.64 jum)/(Br7') » 1 is interpreted as the case of 
shock excitation, as radiative shocks destroy grains and return 
Fe to the gaseous state. In SNR ejecta, however, we do not know 
what fraction of Fe might be in grains initially, or what the Fe/H 
abundance rati o is, making th i s test i nappropriate. In addition, as 
pointed out bv lGraham et al.l (1 19891) . in the Crab Nebula condi- 
tions leading to high values of ([Fe II] 1.64 jum)/(Br'>') can also 
be found in the photoionized filaments. Thus, determining this 
line ratio as the sole factor discriminating between the two sce- 
narios may not be conclusive. Recent stu dies of SNRs w here 
[Fe II] 1. 64 A/m is detecte d, e.g. Gl 1.2-0.3 (iKoo et al.ll2007l) and 
3C 396 ( iLee et al. I I2009I) . show that iron emission can be ex- 
plained by a radiative shock propagating into the gas with a ve- 
locity ~ 100 km s"'. The iron line-emitting regions are charac- 
terised by high electron densities ~ 10* cm"^ and electron 
temperatures of order a few times 10^ K. This can be con- 
trasted with the conditions found for the photoionized case in 
the Crab Nebula whe re n, ^ 200 cm"^ and 4000 < T, < 8000 K 
dGraham et al.lll989l) . Thus in order to try to discriminate be- 
tween the two scenarios in the case of G2 1.5-0. 9, spectroscopic 
observations of the [Fe II] line-emitting regions are needed. 

5. Conclusions 

In this paper we report the first detection in the infrared band 
of emission emanating from the compact nebula region in 
SNR G2 1.5-0. 9 (~2" radius around the pulsar position) through 
imaging and polarimetry. The toroidal structure of the magnetic 
field in the compact nebula is inferred from polarimetric obser- 
vations in the Ks band. The high degree of linear polarisation to- 
gether with the observed swing of the electric field vector across 
the nebula confirms that the observed radiation is synchrotron 
emission from relativistic particles accelerated at the wind termi- 
nation shock and radiating in the downstream flow. The infrared 
spectrum of the compact nebula can be described as a power law 
with spectral index ccir - 0.7 ± 0.3; when combined wi th X-ray 
results from Chandra observations (ICamilo et al.ll2006l) . it sug- 
gests that the spectrum flattens between the infrared and X-ray 
bands. 

The detection of [Fe II] 1.64 jum emitting material within 
the PWN of G2 1.5-0. 9 is also reported. This iron emission is de- 
tected in the NE and in the SW of the PWN. The determined sur- 
face brightness of the SW emission is weaker compared to that 
of the NE region. It is important to point out that the detected 
[Fe II] emission has the shape of a limb-brightened broken ring 
(~40" radius), which follows the edge of the X-ray and radio 
PWN, but is concentrated inside the PWN boundary. Most prob- 
ably the observed iron emission originates from cold supernova 
ejecta shocked by the expanding PWN bubble. In this case col- 
lisionally excited iron radiates in the medium behind the shock. 
However, to confidently discriminate between the case of shock 
excitation and photoionization by the synchrotron continuum, 
spectroscopic observations of the iron line-emitting material are 
needed. Line profiles can give shock velocities, while observa- 
tion of additional spectral lines can provide density and temper- 
ature diagnostics to constrain shock models. Important informa- 
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tion to be gained this way can include the speed of shocks caus- 
ing excitation and ionization; the expansion rate of the PWN; 
abundances of the shocked material; and possible contributions 
of dust to the IR continuum, allowing inferences on newly syn- 
thesized dust. G2 1.5-0. 9 thus joins the small elite class of ob- 
jects, along with G54. 1-1-0.3 and B0540-69.3, where IR observa- 
tions of PWN interactions are illuminating the deep interior of a 
core-collapse supernova. 

The present work demonstrates that infrared observations, 
both imaging and polarimetry, are an effective tool for detecting 
optically obscured pulsar wind nebulae and for studying their 
inner structure. 
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